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Abstract 

The classical Dcbyc-Hiickcl screening effect of the electrostatic field generated by isolated charged 
particles immersed in a plasma is reviewed. The validity of the underlying mathematical model, 
and particularly of the weak-field approximation, are analyzed. It is shown that the presence of 
the plasma sheath around test particles and the resulting effect of charge screening are essential 
for the description of plasmas which are strongly coupled. 
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I. INTRODUCTION 



In this work we intend to analyze the consistency of the traditional mathematical model 
for the so-called Debye screening problem (DSP) originally formulated by Debye and Hiickel 
(DH model In particular, we intend to prove that the 3D Poisson equation involved 
in the DH model does not admit physically acceptable solutions, i.e., solutions which are 
provided by ordinary functions and are at least continuous in the domain of existence, i.e., 
are so-called classical (or strong) solutions. For this purpose a modified model is proposed 
which takes into account the effect of local plasma sheath (i.e., the local domain near test 
particles where the plasma must be considered discrete). Basic consequences of the model 
are discussed, which concern the asymptotic properties of the solutions determined both 
for weakly and strongly-coupled plasmas. As an application the charge screening effect in 
strongly-coupled plasmas is investigated and an explicit expression of the effective charge 
for the asymptotic DH potential is determined. 

An ubiquitous property of plasmas, either quasi-neutral or non-neutral and weakly or 
strongly-coupled, is the so-called Debye shielding of the electrostatic potential. This is 
generated by a single charged particle (test particle) of the plasma, located at position 
r(t) and belonging to the N-body system (with iV ^> 1) of charged particles forming a 
plasma, or more generally by a finite collection of test particles (with positions i"j(t), i = 1, 
N t , possibly with N t <C N). The test particle and the plasma are both assumed non- 
relativistic when referred to a suitable Galileian reference frame, i.e., defined in such a way 
that |V| /c <C 1, r(t) /c€l,V and r(t) being respectively the fluid velocity of the plasma 
and the test particle velocity r(t), both to be assumed generally non- vanishing. It is well- 
known that this phenomenon occurs, in particular, if the plasma is assumed suitably close 
to kinetic equilibrium, namely if each particle species (s) is described by a Maxwellian ki- 
netic distribution function (/ms)? carrying a finite number density N s , temperature T s and 
flow velocity V. Here, T s , V and N s (the latter only in the absence of test particles) must 
result constant or at least, in a suitable sense slowly dependent with respect to position (r) 
and time (t). As a consequence, sufficiently far from each test particle, in a sense to be 
defined, the Coulomb electrostatic potential results effectively screened and decays expo- 
nentially on a characteristic scale specified by the Debye length (X^see below). Of course, 
this phenomenology applies only provided local perturbations of the kinetic distribution 
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function are suitably small. In such a case, due to spherical symmetry of the test particle, 
the screened potential must result manifestly function only of the distance p = |r — r(i)| 
between the point of measurement r and the position of the test particle r(i). In addition, 
if both the test particle and the plasma are non-relativistic, the potential must be also in- 
dependent of the velocity of the test particle r(t). This is the subject of the so-called Debye 
screening problem, which regards the determination of the electrostatic potential generated 
by test particles in a quasi-neutral plasma. The mathematical model of DSP in his original 
formulation jjlj (recalled in detail below, see Sec. 2) concerns point-like test particles and 
is based both on the neglect of the local plasma sheath around each test particle (or the 
finite size of the test particles, for example in the case of large-size dust particles) and on 
the so-called weak-field approximation. The first assumption, concerns the neglect of the 
screening effect produced by the presence of discrete interactions close to the test particles 
(local plasma sheath). This set is generally species-dependent and is defined as the spherical 
domain R{ Pos } = {p : p v < p < p os , P = |r — r(t)\ G 1R + } , p p being the radius of a spherically- 
symmetric test particle and p os the characteristic radius of the plasma sheath produced bt 
the s— th particle species. The characteristic radius p os can manifestly be identified with 

/ „ \l/3 

the mean interparticle distance in the plasma for the species s, i.e., p os = ( AitN J , N os 
denoting the s— th species plasma number density in the absence of test particles. This set 
is centered at the position of the test particle, r(i), in which the species s of the plasma 
must be treated as a discrete system; instead, for any distance p > p os the same species is 
assumed to be described in terms of a continuous kinetic distribution function. The second 
approximation, which permits the linearization of the Poisson equation for the electrostatic 
potential, is based on the requirement that, for arbitrary p = |r — r (t) | G M + , the following 
asymptotic ordering holds 

$(p) ~ 0(e) « 1, (1) 

e being an infinitesimal and $ = |e$| /T e the normalized electrostatic potential and T e the 
electron temperature measured in energy units. As a result, Debye and Hiickel obtained, for 
the electrostatic potential generated by a point-particle carrying the electric charge q, the 
asymptotic solution 

$(p) = * (p) es 47r|r ! r(t)| e-'-W'^ (2) 



3 



which is the so-called DH potential, also known as the Yukawa potential. Here c is the DH 
effective charge to be approximated 

c = g, (3) 

and Xd is the Debye length 



47rZ 2 e 2 iV j 



where the summation is carried out over all plasma species j. Furthermore, A^- — , 
is the j — th species Debye length and N j the j — th species plasma density defined in the 
absence of test particles (i.e., letting $(p) = 0). 

Despite previous attempts to construct approximate or exact solutions to the DH model 
y,Q], the related mathematical model appears incomplete and can be shown to be phys- 
ically unacceptable, due to the neglect of the local plasma sheath. In fact, it is obvious that 
sufficiently close to the point-particle the weak-field approximation (0) is violated making 
the DH model invalid. In the past [3J] it was pointed out that in such a case the test particle 
does not produce any electric field, but only complete charge neutralization by the plasma, 
thus producing a Debye length which effectively vanishes. Other objections concerned the 
asserted indeterminacy of the solution for x = due to its divergence in the same point 
These issues were later addressed in a more general context jj], including the 2D case 
where complete neutralization cannot be achieved. To recover the correct physical picture 
the effect of local plasma sheath must be included. Nevertheless, for suitably dense plasmas 
or in the case of plasma species characterized by very high electric charges (high-Z), such 
as dusty plasmas, the weak-field approximation may be locally violated. This circumstance, 
when the effect of finite local plasma sheath is included, occurs if the normalized electro- 
static potential $(p) results of order unit or larger on the boundary of the plasma sheath 
(produced by at least one of the s plasma species), namely for p = p os . Such a condition 
can be expressed in terms of the Coulomb coupling parameter T^, i.e., by the requirement 
that for a test particle of species k there results at least for a species s 

r _ (3k Z k e 2 

r fcs = — = i (5) 

%os J- epos 

(strong coupling condition). Here the notation is standard. Thus, /3& = ^ , p os = 

( „ \V3 

( 47! /jy ) and x os = Pos/Xd, denote respectively the dimensionless electric charge carried 
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by a test particle of species k and the radius of the s — th species plasma sheath surrounding 
the test particle. We stress that, ordinarily, in plasmas the asymptotic condition x os <C 1 has 
also to be imposed on each particle species s. For example, in a space dusty plasma typical 
values for dust grain charges, plasma temperature and plasma density are ~ 10 4 -f- 10 5 , 
T ~ 1 eV, N ~ 10 9 cm -3 . In this case the Coulomb parameter for a negatively-charged dust 
particle in the presence of the plasma sheath produced by hydrogen ions may result typically 

Z ( — N cm -3 )'' 3 

T di = 0.91 x 10~ 6 n 3 (r eV ) — — = 47r x (15 + 150) , while the dimensionless radius of the ion 

(No cm" 3 ) 1 ' 6 

plasma sheath can be estimated x Q i = 0.83 x 10~ 3 v - ' ^-^j — — 0.03. However, in principle, 
in a dusty plasma the plasma Debye length can also become significantly smaller than in the 
corresponding "clean" plasma (i.e., before injection of the dusty species), which may influ- 
ence, i.e., increasing, the values of and x {. Further important aspects of DSP concern the 
investigation of possible charge screening effects, namely the reduction to the effective charge 
c with respect to the asymptotic value characterizing the weak-field approximation Q of 
the DH potential (0). These effects, produced by non-linear contributions in the Poisson 
equation, while usually negligible for weakly coupled plasmas, are known to be significant in 
strongly-coupled plasmas pj], such as dusty plasmas U, 3|. Indeed the investigation of the 



effective interactions characterizing high-Z grains in plasmas has attracted interest in recent 
years especially for numerical simulations (see, for example, 12, 14, ^|)- However, analytic 



estimates of the effective charge characterizing the DH potential in strongly-coupled plasmas 
are still not available. 

Goal of this work is the analysis of DSP and the definition of a suitably modified math- 
ematical model to take into account the effect of local plasma sheaths in quasi-neutral 
plasmas. In particular, in Sec. 2 the traditional DH model is recalled, while the new model 
is presented in Sec. 3. In the same section the basic mathematical results are given which 
concern the asymptotic properties of the solutions obtained in the limit of vanishing sheath 
radii (p os — > + ). We intend to prove that the Poisson equation characterizing the DH 
model (here denoted Debye-Poisson equation) can be considered as a limit equation ob- 
tained from a suitably modified formulation of DSP. The latter is obtained by including the 
effect of finite local plasma sheaths (i.e., requiring x os > 0, for all species). Basic feature of 
the present approach is the representation of the Poisson equation in integral form. This 
permits to analyze the asymptotic properties of the solutions of the modified problem in the 
limit x os — > + . It is found, that the limit solution of the modified DSP for x os — ► + is a 



distribution which vanishes identically for all p > and is discontinuous in p = 0. It is found 
(see THM.l in Sec. 3) that in this limit the solution of the modified DSP results uniquely 
determined as a function of the normalized charge carried by the test particle (/?) and is 
represented by a well-defined distribution. Precisely it follows that (see THM.l in Sec. 3) 

lim smh$ Xo (x) = -^S(x), (6) 

x o ^0+ X z 

while one can prove that the limit function lim Xo ^ + &x ( x ) = 3>(#) reads (see THM.3 in 
Sec.4) 

lim $ Xo (x) = lim $ Xo (x ) exp {y {int) (x, x Q ) } , (7) 

x- o ^0+ x Q — +0+ 

where y^ nt ^(x, x Q ) is a suitable smooth real function to be defined [see Sec.4, Eq. (ffiUj) ] so that 
there results y(x ,x ) = while for x ^ x Q the following limit holds 

lim y iint \x,x ) = -oo. (8) 

x o ^0+ 

Hence, obviously, the limit function <3>(x) is not a strong solution of the DP equation. This 
is therefore a characteristic property of the DH model. In particular, as a basic consequence, 
the effective charge of the DH asymptotic solution c vanishes identically in such a limit and 
results independent of the charge of the test particle. 

As a further development, we intend to investigate strongly-coupled plasmas (Sec.4 and 
5), for which there results for a test particle of species k (for example to be identified with 
a dusty grain) 

r * ~ W) >> 1 (9) 

x c ,~0(S t ), (10) 

(strong coupling ordering) where k — 0,1 and 5 is an infinitesimal. The asymptotic solutions 
of the modified DSP are here determined explicitly for the external and internal asymptotic 
solutions. The former, which are valid outside the Debye sphere (i.e., in the domain x > 1, 
or x ^> 1) coincides with the customary DH potential. The latter instead, which occurs 
inside the Debye sphere near the boundary of the local plasma sheath, describes the non- 
linear effects due to the screening generated by plasma species having opposite charge with 
respect to that of the test particle. In addition, in Sec.4 also the asymptotic solution of 
the modified DSP close to the boundary of the local plasma sheath (i.e., for Ax = x — x os 



infinitesimal) is constructed. It is proven that the asymptotic solution can be determined 
uniquely, together with its internal boundary conditions (defined on the boundary of the 
plasma sheath, i.e., for x = x os ). Finally, by comparing the internal and external asymptotic 
solutions for a test particle in a strongly-coupled plasma, an "asymptotic" upper bound is 
obtained (Sec. 5) for the effective charge characterizing the DH potential. As a result, it is 
found that in strongly-coupled plasmas the effective charge of highly-charged test particles, 
which characterizes the DH asymptotic solution, results strongly reduced with respect the 
value of the isolated test particle (charge screening effect). The result appears potentially 
relevant especially for strongly-coupled dusty plasmas, since the charge screening effect can 
dramatically reduce the DH potential, particularly for highly charged test particles. 



II. THE DEBYE-HUCKEL MODEL AND THE DSP PROBLEM 

The traditional formulation of the DSP, based on the Debye-Hiickel model Q regards 
the test particles as point-like and having a spherically-symmetric charge distribution while 
ignoring the effect of local plasma sheath. This implies, from the physical standpoint, to 
neglect the discrete nature of the plasma. As a consequence the background plasma is 
treated as a continuum medium formed potentially by several particle species, each de- 
scribed by a Maxwellian kinetic distribution, so that its total charge density results as 
p p = ^2 Z s eN os exp < — 4pf^ r > where the constraint Z s eN os = must be imposed for 
quasi-neutral plasma. Here the summation is performed on all particle species of the plasma, 
while N os and T s (species density coefficient and temperature) are both assumed constant. 
Furthermore, the charge density of the test particles is manifestly p n = ^2 rrnlr ) where 

47r|r— rj(t)| 

q,i = Z q .e and r^(t) are the electric charge and position vector of the i— th test particle and 
the summation is carried out on the collection of all test particles. In the remainder we shall 
assume, without loss of generality, for all test charges qi > 0. Hence, in terms of dimen- 
sionless electrostatic potential <3> = ^j?, one obtains, for example, p p = 2eA sinh$, for a 
quasi-neutral electron-proton plasma. In a case of a single point charge, in terms of the De- 
bye length \ D = and the dimensionless distance x = p/Xd (where p = |r — r(t)| g), 
$(x) must satisfy the equation: 

V^$(x) = + sinh8(x), (11) 
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(DP equation), where V 2 , = -^-^ (x 2 -^) (the Laplacian expressed in dimensionless vari- 
ables). We require that <&(x) obeys the boundary conditions, i.e., 

lim 8(x) = 0, (12) 



x— >oo 



\imx 2 -^-$(x) = -13. (13) 
x^o dx 

Equation ([lip, together with (|12p. ()13|). define the so-called Debye screening problem (DSP). 
We notice that, in contrast to what stated earlier p], both boundary conditions ()12)) . |T3|) 
can be satisfied by $(x). In fact, it is immediate to prove that in the domain x > 0, and 
imposing ()12j) . ()13j) . the DP equation can be cast in the integral form 



$(x)=- 

X 



^ rx poo 

— / dx'x' 2 + / cfccV 
x Jo 



sinh$(x / ). (14) 



A strong solution of (|TT|) . (fT2^ . (JTHJ) [or, equivalent, of Eq.((T3J)] is a solution belonging to the 
functional class <5 (2) (M + ) = : C (0) ([0, oo[) , C (2) (]0, oo[)|. In the traditional approach 

the DSP problem is solved requiring the validity for all xG + of the weak electric field 
ordering ((TJ) and assuming that $(x) =& (x) + 0(e) (weak field approximation), § (x) being 
the solution of the linearized DP equation : 

V 2 3 (x) = -^5(x) + ^ (x), (15) 



x 2 

with & (x), by assumption, belonging to the same functional class of Q(x). This equation 
admits the exact solution § (x) = fe~ x , i.e., Eq.(0), defined in terms of the DH effective 
charge c = j3. Finally, we stress that, assuming the validity of the weak electric field ordering, 
due to the linearity of Eq. (jl5)l a similar result is implied also for a collection of point charge 
sources. In this case, in fact, the corresponding solution of (|TK|) is just a sum of DH potentials, 
one associated with each source. 



III. CONSISTENCY OF THE MATHEMATICAL MODEL - THE MODIFIED DSP 

In this section we define the modified Debye screening problem, based on the introduction 
of a simple model of local plasma sheath. In the sequel we shall consider for simplicity of 
notation the case of a two species-plasma, formed by electrons and Hydrogen ions, having an 
unique plasma sheath. However, the generalization to multispecies plasma is straightforward. 
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Thus, we shall assume that the test particle is represented by a spherically symmetric charge 
of radius p p . For a particle in which p p < p a the plasma sheath is represented by the spherical 
shell centered at the position (center) of the test particle for which p p < p < p , in which 
the plasma charge density (except for the presence of the test particle) results negligible. In 
the sequel we can also let in particular p p = (point-like test particle) or p p = p Q (finite-size 
test particle). 

The customary DH model is obviously recovered letting p p = and taking the limit 
p — > (or in dimensionless variables, requiring x p = p p /\o = and x Q = p /^D —> 
0). Denoting § Xo (x) the solution of the Poisson equation, here we intend to determine its 
asymptotic properties in the limit Xq—* + , while also letting x p = (see THM.l). As a 
consequence and in agreement with 0, Q , in such a case it follows that the limit function 
\im Xo ^Q+ Q Xo (x) = $(x) vanishes identically for x > 

lim $ Xo (x) = 0. (16) 

In addition, in the same set we intend to prove the identity 

(3- lim / dx'x' 2 smh$ Xo (x')Q(x' - x Q ) = 0, (17) 

*o^0+ J Xo 

which implies necessarily Eq.®, Q(x — x a ) being the weak Heaviside function 

e(x-x ) = { (is) 

I 1 X > x _ 

In detail the relevant equations valid in each subdomain for the normalized electrostatic 
potential $ Xo (x) are as follows. In the internal domain < x < x p the electrostatic potential 
is assumed constant 

®xM =$ Xo (x p ). (19) 

In the plasma sheath x p < x < x p , § Xo (x) satisfies the customary Poisson equation in the 
presence of the charge density produced by a finite-size spherically-symmetric charge 

Vl$ Xo = -^S(x-x p ). (20) 

x A 

Finally, in the external domain x > x there holds the Poisson equation in the presence of 
the plasma charge density: 

V 2 x $ Xo = e(x-x )smh$ Xo . (21) 



The boundary conditions analogous to (jU j) . (fT5 j) . imposed respectively at infinity and at the 
boundary of the plasma sheath, are specified as follows 



x—*oo 



lim$, o (x) = 0, (22) 
= -P. (23) 



We notice that, if x p < x a (for example, x p = 0), & Xo (x) results by assumption at least 
of class CW(M{ Xp }), where ^{ Xp } = ]x p ,oo[ . Here x , (3 are both assumed constant and 
strictly positive real numbers. The problem defined by (|2"T|) . (f20|) . together with the boundary 
conditions (|2*2*)l . (|2*3|) . will be here denoted as modified DSP. From the physical standpoint 
Eqs. f{21j) .(pn jl may be viewed as the Poisson equation for a spherical ideally conducting 
charge, or for a point particle in the presence of a plasma sheath, of radius r Q (i.e., x Q = r Q /\ D 
in non-dimensional variables) which is in electrostatic equilibrium and is immersed in a 
spatially uniform quasi-neutral and Maxwellian plasma. As for the previous DP equation, 
it follows that, for solutions satisfying the boundary conditions (J22|) . (|2"3*j) . in the domain 
x G Rra; } Eq. (j21j) can be cast in the integral form 

$iAx)J m^A_ (24) 

x 

' FX POO 

— / dx'x' 2 + / dx'x' sinh <& Xo (x')0(x' — x a ). (25) 

. J x J X 

In particular, thanks to continuity at x = x Q of $ x . o (a;), one obtains the constraint 

/•oo 

$ Xo (x ) = T- / dx'x' sinh $ Xo {x')Q{x' - x ), (26) 

J X 

with r = ^- denoting the Coulomb coupling parameter. It is immediate to establish, the 
existence and uniqueness of & Xo (x) in the functional class C^°°' ) (lR{ 2:o }), together with its con- 
tinuous dependence on initial data, in particular the continuity with respect to the parameter 
x a G{o} • Moreover, assuming that the weak- fields approximation (JTJ) applies (this condition 
is manifestly fulfilled identically in the weak-coupling ordering, < T ~ 0(e) 1), and is 
satisfied at least for x 3> 1 suitably large, it is immediate to prove that in this subset an 
asymptotic solution of the modified DSP is provided by the external asymptotic solution 



$ Xo (x) = ® x e o xt) (x) = -e~ x+x °. (27) 
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Here denoted as of the modified DSP and c = c(x a , T) is the effective dimensionless charge. 
Hence, $> Xo (x) reduces formally to the previous DH potential (j2J) when x ^> x Q . In the 
weak-coupling ordering it follows c(x 01 T) = -rf— while for strongly-coupled plasmas a 
lower value is expected. Furthermore, it is obvious that the limit function lim Xo ^o+ &x a (%) 
coincides with the solution of DSP, i.e., 

lim $> Xo ( x ) = ®( x )- ( 28 ) 

Let us now investigate the asymptotic properties of the exact solution $ Xo (x) in the limit 
x -> 0+ 



A. THEOREM 1 - Asymptotic properties of $ Xo {x) 

For any strong solution of the modified DSP, & Xo (x) obtained letting x p = 0, the limit 
function \im Xo _ + $ Xo (x) has the following properties: 

lim $ Xo (x ) = +oo; (29) 

Xo— >0+ 

2) for any x > Eg. Mb}) holds identically. 

3) for any x > the integral limit \11\J is satisfied by & Xo (x).This implies that the limit 
function = \im Xo ^ + Q Xo (x) results such that Eq.fiQ) is satisfied identically. Moreover, 
for any x > 0, x G{o} 

$0) = lim § Xo (x) = 0. (30) 

x o ^0+ 

4 ) the following limit is satisfied by the boundary value & Xo (x Q ) 

lim x $ Xo (x ) = 0. (31) 

x o -*0+ 

5) the limit value of the effective dimensionless charge c(x Q , T) for x Q — > + , obtained 
keeping T finite, is 

lim c(x ,T) = 0. (32) 

x o ^0+ 

PROOF 

1) In fact, as a consequence of the integral equation (}2l)j) and the continuous dependence 

of $ Xo (x) on the initial data, it follows 

^ f°° ^ ^ 

lim x Q Xo (x ) = (3 — lim x Q I dx'x' sinh <& Xo (x')@(x' — x a ), (33) 

x o ^0+ x o ^0+ J x 
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which implies 



lim $ Xo (x ) = oo, (34) 

x o ^0+ 

POO 

lim / dx'x' sinh Q Xo (x)Q(x' — x Q ) = oo, (35) 

x o ^0+ J Xo 

i.e., the limit function lim x . o ^ + &x a ( x ) diverges in x = x Q . Therefore, due to the continuity 
of §x (x) with respect to x G [0, oof it follows that infinitesimally close to x,x Q = 0, and 
when infinitesimal of the same order, Q Xo (x) must diverge logarithmically as 

KW-ln^y (36) 

2) Let us now consider the implications of the integral equation (|24|) for the limit function 
lim Xo ^o+ ^x (x) for arbitrary x G{x } • There follows 



^ 3 1 f x - - 

lim ® x (x) = — lim / dx'x' 2 sinh $ x (x')6(x' - x )- (37) 

poo ^ ^ 

— / dx'x' lim sinh $ 1 . o (x / )0(x / — x Q ), 

Jx x o ->0+ 

where, due to the asymptotic estimate (j3BJ), the second term on the r.h.s. necessarily diverges 

1 r x ^ ^ 

lim — / dx'x' 2 sinh $ Xo (x')Q(x' — x a ) = 00 (38) 

unless there results for any x 7^ x , x &{ Xo } 

lim $ Xo (x) = 0. (39) 

3) As a consequence of Eq. ()39j) . from the integral equation ()24)1 it follows necessarily that 
for all x > 0: 

lim ^— — = lim — / dx'x' 2 sinh <& x (x')O(x' — x Q ), (40) 

x ^o+ x s ^o+ x y a . o 

which proves the limit ([17)1 . As a consequence it must result necessarily that the limit 
lim :ro ^o+ sinh$ Xo (x) is a Dirac delta, in accordance with Eqs. (|16|) and (JHJ). 

4) The limit (j31|) follows immediately from the boundary condition (|26jl. 

5) Finally, Eq. ()39|l implies manifestly the limit ()32|1 . 

As an immediate consequence of THM. 1 it follows that the DP equation of the DH model 
does not admit strong solutions. 

12 



B. THEOREM 2 - Non-existence of classical solutions of DSP 



In the functional class C^ 2 ^({o}) the DSP problem defined by Eqs.\ll )) .\l ty) and A ISA) has 
no strong solution. 



PROOF 



In fact, first, we notice that the limit function 



lim $ Xo (x) = $(x), 



(41) 



is manifestly a solution of the DP equation which satisfies the required boundary conditions 
(JT2J, (fT3|). Other hand, due to THM.l, this solution is discontinuous in x = x and results 
a distribution. Hence it is not a strong (classical) solution of the modified DSP problem. 

The basic implication of THM.l and 2 is that the DP equation ([lip, provided by the 
DH model, must be regarded as physically unacceptable, since it does not admit strong 
solutions. In this regard it should be noted that, as a basic principle, physically acceptable 
of solutions of ordinary (or partial) differential equations characterizing the classical theory 
of fields must be suitably smooth strong solutions. The modified Debye screening problem 
here defined, instead, exhibits smooth strong solutions and therefore appears, from this 
viewpoint, consistent. 

IV. INTERNAL ASYMPTOTIC SOLUTION OF THE MODIFIED DSP 

In this section we intend to construct explicitly the internal asymptotic solution of the 
modified DSP problem, defined by (|2i p .(|2D jl . (J22J) ? (E3) , in the case of strongly-coupled plas- 
mas, for which the asymptotic ordering defined by the inequalities (J9j). (fTUJl applies. In 
particular, we shall consider a suitable neighborhood of the boundary of the plasma sheath 
(x ), characterized by an infinitesimal amplitude Ax = x — x Q defined so that 



where by hypothesis w Q = In y^rj ^ 1- Under such assumptions, we intend to construct 
approximate solutions of the modified DSP which are accurate to leading order in the in- 
finitesimal dimensionless parameter 5 and with respect to Ax. These will be denoted as 
internal asymptotic solution of the modified DSP. We intend to prove, as a consequence, 



< Ax < 




(42) 
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that in the neighborhood of x Q defined by (|42|) the internal asymptotic solution can be 
uniquely determined, to leading order in the relevant asymptotic parameter. In fact, the 
following theorem holds. 

A. THEOREM 3 - Internal asymptotic solution for strongly-coupled plasma 

In validity of the asymptotic ordering of strongly- coupled plasma [see inequalities Hty. AlCty J. 
requiring x p = (or at least x p < x D ) and in the neighborhood of x Q defined by the 
inequalities the asymptotic solution of the modified DSP [$EM),$£W> reads 

$ Xo (x) = *W(x) = $ exp{y( int \x,x )} (43) 

where $„ = $> Xo (x ) and 

y {int) (x, x ) = y (x', x ) = -a(x- x Q ) (44) 

(internal asymptotic solution), where a = ^i- and the integration constant $ Q results 

/2F 2 \ 

•• sb, UJ- (45) 

PROOF 

In order to prove these results let us represent $> Xo {x) in the form 

§ Xo {x) = $ exp{y{x,x )} , (46) 

It follows that in the domain x > x Q , y(x,x D ) satisfies the ordinary differential equation 

2 1 r ^ i 

y" + {y'f + ~y' = exp <^ $ D exp {y} - y(x, x ) \ - (47) 
x 2$ ^ J 

^exp {-8 exp{y} -y(x,x )\ = S(y(x,x )), (48) 

2$ c L 

where primes denote differentiation with respect to x. Consistent with (|22 |) .([23 )) . there follow 
the boundary conditions 

y(x o ,x o )=0, (49) 



lim y(x, x ) = — oo, (50) 
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lim y(x, x Q ) = — oo, 

x— >oo 



(51) 



v'M\ x=Xo = (52) 

As a consequence, it is immediate to prove that y(x,x ) for arbitrary x,x Q > (x > x D ) 
with admits the exact integral representation 

y{x,x ) = y (x,x ) + F(x,x ), (53) 
y {x', x ) = -a(x - x ). (54) 

Here F(x, x Q ) is the solution of the integral equation 

X 

F(x,x a ) = J dx' (x - x') ^S(y(x',x )) - y(x',x ) 2 - -^y' (x' , x )j , (55) 

Xq 

while S(y(x', x )) is defined by the r.h.s. of Eq. (}4Tj) . Furthermore, a direct evaluation of the 
explicit integral in Eq.()55|) permits to cast F(x,x ) in the form 

F(x,x ) = F (x,x ) + Fi(x,x ), (56) 

X 

^) = S ix ' (x - x ,)GM (57) 

X o 

where F (x,x ) and Gi(x',x Q ) are respectively 

X 

F (x,x ) = J dx' (x - x 1 ) ^S(y (x',x )) - y (x',x ) 2 - -^y (x',x )\ , (58) 

Xo 

2 

Gt(x' \x ) = S(y(x',x )) - S(y {x',x )) - Fi(x',x ) 2 - 2F 1 (x',x )y (x',x ) -F{(x',x ). 

x 

(59) 

Using Eqs. (J3HI), (fSHj) and (JBTJj) it is immediate to prove that y(x,x ) satisfies also the 
boundary conditions (}4T?|) . (JBTj) and (|32*j) and in particular that the limit (JSTij) is actually 
satisfied. The proof follows by noting that, in the strong-coupling ordering (}9*|).()10J) and for 

< (x — x ) ~ 0(5°), the dominant contribution in y(x,x ) is negative and is provided by 
the second term on the r.h.s. of Eq. (J58j) . This implies that an asymptotic approximation 
for y(x, x ) is provided in this case by: 

y(x, x ) = y {int \x, x ) = -a(x - x ) + F a (x, x ), (60) 
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and consequently also Eg . (J44j) . Manifestly y( mt >(x,x ) satisfies identically all the boundary 
conditions (jUIJ), (JH0*|) . (J5T|) and (JH2J) . Finally, to determine the initial condition § Xo (xo) = $0, 
it is sufficient to notice that in asymptotic orderings (|9|).([TU]) and (j4"2"j) the following identity 
must be satisfied, correct to leading order in 5, 

P~\j x dx'x' 2 exp{$ exp{y(x',x )}}. (61) 

The integral can be estimated asymptotically to yield 



- 2T 
$x (x ) = In 



r-^ o (x ) 



(62) 



which to leading order in 5 implies Eq.()45|). To briefly comment these results, it is important 
to remark, that, by assumption, the internal asymptotic solution here obtained, ([60jh is 
valid if the parameter x Q results of order 0(1) or is infinitesimal in S [see the ordering ([10))]. 
while the parameter T satisfies Q. As a consequence, in the limit x a — > + Eq. fJBTj) results 
consistent with Eq.(J7|), thus proving also the validity of Eq.fjHJ). Finally, we notice that the 
sub-domain (|42j) in which the internal asymptotic solution (|6U|) is valid, allows & Xo (x) <C 1. 
This is also consistent with the weak-field approximation (JTJ). 



V. CONCLUSIONS: ASYMPTOTIC UPPER BOUND FOR THE EFFECTIVE 
CHARGE 

As an important application, the previous analytic results can be used to obtain an 
asymptotic estimate for the effective dimensionless charge c(x Q , T) carried by the DH poten- 
tial in strongly-coupled plasmas. We intend to show that test particles which have a high 
charge undergo a strong charge screening effect, close to the local plasma sheath, produced 
by non-linear effects in the Poisson equation. As a consequence, outside the Debye sphere 
(i.e., for x > 1) the DH potential generated by highly charged test particles in strongly- 
coupled plasmas results strongly reduced with respect to the theoretical value observed in 
the corresponding weakly-coupled systems. 

The effective charge of the test particle can be estimated by comparing the internal and 
external asymptotic solutions in a suitable neighborhood of the boundary of the plasma 
sheath (i.e., near x = x Q ), belonging to the sub-domain defined by (j4"2"j) . For this purpose 
we impose the matching condition 

$r) ( * c )=$rw, (63) 
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where, by assumption, xc is defined so that there results Q Xo (x c ) <C 1 and moreover 

' Ax c ~ 0(5) < 1. (64) 



being Ax c = x c — x Q and where — ~- ^> 1, consistent with the assumption of strongly- 
coupled plasma ([§ ]) . (|10 p . The matching condition can be used to obtain an upper estimate 
for c(x G ,r). In fact, in validity of the ordering Eq. (j64"|) . $i„ (x c ) can be approximated 

by $a"^(x c ) = $ exp < — x Q ) \ . Then, ignoring higher order corrections in 5, the 

following inequality manifestly holds: 

c(x ,T) < cW(x ,T) = x \n—, (65) 

x 2 

which delivers therefore an asymptotic upper bound for the effective charge c(x ,T). In 
particular, it follows that T) depends logarithmically on the normalized charge of 

the isolated test particle (3. Thus, in general for strongly- coupled plasmas the effective 
dimensionless charge appears much smaller than in weakly- coupled plasmas. The result can 
be significant for the investigation of dusty plasmas, particularly to describe the charge 
screening effect of highly charged dusty grains. Comparisons with numerical simulations, 
based on the numerical solution of modified DSP [defined by Eqs. (|2T j) .(pD |) and the boundary 
conditions (J22 }) .p5 j) ]. indicate that the asymptotic estimate c^ a ^(x ,r) holds already for 
/3 > 1 and x Q < 0.3. For numerical calculations, the effective charge c(x a , T) has been 
determined numerically as the limit for x ^> 1 (i.e., x — > oo) of the the function 

Ux) = (66) 



where $ Xo (x) denotes the numerical solution of the modified DSP. As previously discovered 
[5], it is found that for strongly-coupled plasmas (i.e., assuming T > 1) f c (x) approaches 
rapidly the asymptotic value c(x Q , T), even at a "distance" x smaller than the Debye Length 
(i.e., for x < 1). For reference in figs 1-3 the numerical estimates of f c (x),c(x ,T) and 
c^ a ^(x ,r) have been reported for the cases j3 = 1,5, 10 and x Q = 0.05 -r- 0.3. Its is found 
that the upper bound provided by the majorization (}65|) can actually be used to 

obtain an approximate estimate of the effective charge c(x ,T). 

Finally, it is worthwhile pointing out that the asymptotic estimate provided by c^ a ^(x , T) 
satisfies also the correct limit set by THM.l. In fact the it results 

lim c(x , T) = (67) 
17 



and this limit is satisfied even in the case in which T ~ 1/5 ~ l/x a . 
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Figure captions 



Figure 1 - Comparison between f c (x), c(x ,T) and c^(x 0: T). The data are normalized 
with respect to (3, the normalized charge of the isolated test particle. The figure concerns 
the case with (3 = 1 and x Q = 0.05, yielding T = 20. The horizontal straight line represents 
the asymptotic estimate c < - a ^(x D , T), while the curve below it is the graph of f c {x). It follows 
f c (x ) = 0.564, while the asymptotic value c(x Q , T) = 0.493 is reached approximately at 
x ~ 0.3, and the upper bound for the normalized effective charge is c^(x Q , Y) = 0.589. 

Figure 2 - Comparison between f c (x), c(x Q , T) and c^(x Q , T) for (3 = 5 and x Q = 0.2 (with 
T = 25). In this case f c (x ) = 0.369, while the asymptotic value c(x ,T) = 0.28 is reached 
approximately at x w 0.4, and c (a)( (a; ,r) = 0.38. 

Figure 3 - Comparison between f c (x),c(x ,T) and c (a)( (a; ,r) for (3 = 10 and x Q = 0.3 
(with T = 33). In this case it is found f c (x ) — 0.273, while the asymptotic value c(x ,T) 
= 0.188 is reached approximately at x m 0.5, and c (a)( (x ,r) = 0.303. 
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